Finite size corrections to the blackbody radiation laws 
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We investigate the radiation of a blackbody in a cavity of finite size. For a given geometry, 
we use semiclassical techniques to obtain explicit expressions of the modified Planck's and Stefan- 
Boltzmann's blackbody radiation laws as a function of the size and shape of the cavity. We determine 
the range of parameters (temperature, size and shape of the cavity) for which these effects are 
accessible to experimental verification. Finally we discuss potential applications of our findings in 
the physics of the cosmic microwave background and sonoluminescence. 
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To a good approximation a blackbody can be repre- 
sented by a cavity heated to a temperature T and con- 
nected to the outside by a small aperture. A useful quan- 
tity to describe the radiation properties of a blackbody 
is the density of energy per unit of frequency, 



(1) 



where V is the volume of the blackbody cavity, p{v)dv 
is the number of stationary electromagnetic modes with 
frequencies between v and v + dv and e{v) is the average 
energy per mode. In the V oo limit. 



p{v) = pv{v) = V 



where c is the speed of light. Classically the average 
energy per mode does not depend on the frequency v 
so, uiy) oc . However experimental results indicate 
uiy) — » 0, — > CXI. Planck noted that agreement with 
experiments could be achieved if the energy was consid- 
ered a discrete variable, 



t{v) = hvj [e^BT - 1^ 



(hs and h are the Boltzmann and Planck constants re- 
spectively). Combining these two results, Planck's radi- 
ation law follows immediately, uy{y) = ^^y*^^ e'"'/'^B^-i ' 
Another useful quantity to characterize the blackbody 
radiation is the total energy emitted per unit of time and 
area, 



R{T) = — y d.vu{i 



In the F — > oo limit it is given by the Stefan-Boltzmann 
law, 



R{T) = 

where a = 5.67 x lO^^Js^^i 



''K-^. The blackbody ra- 



diation, besides its role as a precursor of the quantum 
theory, has been an invaluable tool in a variety of appli- 
cations. For example, it opened the way to the thermal 



remote sensing of bodies that allows to diagnose the prop- 
erties of materials or to measure the surface temperature 
of celestial objects. The international scale of tempera- 
ture, in the range above the freezing point of Silver, is 
based on Eq. 

Finite size corrections to these results are expected 
since pv (j^) is only the leading term of a full expansion in 
powers of a typical length L of the cavity. A comprehen- 
sive theoretical study of p(j^) was carried out in P, Q by 
using the Green function formalism previously developed 
for scalar waves in Refs.0, jij. In the special case of 3d 
cavities with smooth surfaces and L 3> c/v, 

P{v) =y^6{v - v{) ^ pv{v) + PC+ p{v)-, (2) 



where Vi stands for the ith natural mode of the cavity, 
PC = ^ ^ and C oc L is the mean curvature (C = Atvtq for 
a sphere Q of radius tq), namely, the mean of the two 
principal radii of curvature of the cavity [5] integrated 
over solid angle. This oscillatory part is given by, 

p{^') Ap{i', L) cosjvlp/ c + pp) 



where the sum runs over the periodic orbits, p, of length 
L of the classical counterpart. For electromagnetic waves 
P, 01 periodic orbits are the trajectories of the light rays 
inside the cavity as dictated by geometrical optics. The 
amplitude Ap{v, L) and the Maslov index pp can be eval- 
uated explicitly from the knowledge of the classical pe- 
riodic orbits (see Appendix in Ref. [5|). In 3d cavities 



Ap{iy,L) cx (L/c)i+"/2j,"/2 



where n is the number of axes of symmetry of the cavity 
{n — 3 for a sphere). Polarization effects also modify 
this amplitude To leading order only periodic orbits 
with an even number of reflections on the boundary of 
the cavity contribute to the amplitude. Moreover the 
amplitude picks up an additional factor 2 (2cosa) for 
planar (non planar) orbits due to polarization effects. For 



2 



a given periodic orbit characterized by m reflections, 

2 cos a = w ■ w„i + w' ■ w'^ 

where w, w' are two orthogonal unit vectors transverse to 
the direction of propagation of the light ray transforming 
into Wm,w'^a after m reflections. The angle a is usually 
referred to as the polarization angle. 

We note that in this approach p(y) is effectively written 
as an expansion in power of the parameter ^ with A = 
^. The limits of applicability of this formalism are thus 
restricted to a range of sizes L and frequencies v such 
that X ^ 1- 

The main goals of this paper are: a) to provide a de- 
tailed account of the impact of finite size effects on the 
blackbody radiation, b) to determine the range of param- 
eters for which such corrections are accessible to exper- 
imental verification, c) to propose potential applications 
of these results. 



CORRECTIONS TO PLANCK'S AND 
STEFAN-BOLTZMANN'S LAWS. 

Finite deviations from the Planck's law Eq.Q are only 
due to PC and p. As a consequence, the modified Planck's 
law can be written as, 



u{v)_ 

Uv 



1 



„2-n/2 



ai 



L2zy2 (j,L)2-"/2 



Ap cos{lpiy/c + Mp)-(3) 



where ai oc £ is a dimensionless coefficient that depends 
on the mean curvature C oc L of the cavity, Ap = Apc/ L 
is the dimensionless amplitude corresponding to a single 
periodic orbit of length Ip. In the case of symmetric cav- 
ities like spheres (n — 3) and rectangles (n = 2) these 
coefficients can be obtained easily (see below and Q). 
For chaotic cavities (n = 0) it is much harder to get ex- 
plicit analytical expressions for Ap except in the case of 
the shortest periodic orbits 

In order to fully determine we have yet to set a 
cutoff in the sum above. From a physical point of view 
it is evident that periodic orbits longer than the length 
lesc associated to the typical time that a ray stays in the 
blackbody cavity before escaping through the aperture 
cannot contribute significantly to u^iy). The explicit ex- 
pression of the cutoff function depends strongly on the 
symmetries of the cavity. For chaotic cavities {n — 0), 
the contribution to Eq.Q of a light ray of length Ip 
is weighted by the probability that this ray does not 
escape through the aperture, Pesc{lp) ~ 6"'"/'=" with 
iesc ~ ^ Imin, where V is the volume of the black- 
body, Imin is the length of the shortest periodic orbit 
and Aape is the surface of the aperture [1]. By contrast 
for symmetric cavities (n = 2,3) the cut-off function has 
power-law tails Q . 



The Stefan-Boltzmann's R{T) = ctT"' law is also mod- 
ified by finite size effects. After integrating 



R{T) = RiT) + Rc{T) + cjT^ 



(4) 



where 



Rc{T) 



1 



727r L'^h 



[Tkef oc l/L^ 



comes from the curvature term pc in the spectral den- 
sity. The contribution coming from the fluctuating part 
of the spectral density R{T) can be written exactly in 
terms of polygamma functions. This expression is rather 
cumbersome but in the limit Imin > (relevant to 

experiments) simplifies considerably, 



(5) 



where the dimensionless coefficients bo and bi depend on 
the number, n of symmetries axes of the cavity. The case 
n = 3, the sphere, will be discussed in detail later on. For 
n = O.n = 2, 



bo^-- 



2 (-2t:L\ 



l+n/2 



Ap,r{L)f{pp) 
(Tpr)2+»/2 



with f{pp) = cos(/ip) 4- (sin(^p) - cos(^p))n/2. 



bi = 



/2ttL\ 



\ "/2 
' p,r 



Ap^r{L)g{pp) 

(^^^)l+n/2 



with g{pp) = — sin(/ip) — (sin(/ip) — cos{pp))n/2. Sim- 
ilar expressions are found for n = 1. The amplitude 
Ap^r describes the contribution of a single periodic or- 
bit of period Tp = Ip/c repeated r times. We note that: 
a) this sum is convergent even without including a cut- 
off function Pesc{lp) due to the aperture, b) the largest 
contribution to the sum comes from the shortest peri- 
odic orbits (light rays), c) the Maslov indexes pp and 
amplitudes Ap^r for the shortest periodic orbits can be 
computed analytically [S^, d) for T ^ and fixed L our 
semiclassical formalism fails since the maximum of u{v) 
is in the region X/L > 1. 



FLUCTUATIONS IN CAVITIES WITH NO 
SYMMETRY AXES. 

An explicit analytical determination of R(T) or the os- 
cillating part of u(i^) involves the knowledge of periodic 
orbits (light rays) of different lengths. This determina- 
tion is only straightforward in highly symmetric cavities 
with several symmetry axis. However semiclassical tech- 
niques based on that fact that the classical dynamics is 
ergodic [6| permit an estimation of these deviations in 
the case of chaotic cavities with no symmetries. We first 
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study cx Ep,p',.,.'A 



A , , 

■p,r^p' ,r' 



where the average is 



over cavities characterized by the same typical length L 
and with no axes of symmetry. To leading order, the 
double sum Q above is given by the diagonal (p = p') 
term, 



of Fresnel integral that smoothly modulates the ampli- 
tude and phase of the fluctuations. The effect of small, 
non overlapping bumps Ref.[l8| is only to suppress pe- 
riodic orbits of the sphere longer than a certain cut-off 
related to the typical size of the bump. 



4L6 



dr 



K{t) 



L6 



(0) 



where K{t) — J2p r^iT — rTp). The lower limit of inte- 
gration corresponds to the period of the shortest periodic 
orbit Tmin = Imin/c. In the region around Tmin the func- 
tion K{t) is a collection of well separated peaks located 
whose exact positions positions depend on the form of 
the cavity 6]. However for longer periods the number of 
orbits (light rays) with similar periods increases dramat- 
ically (in cavities with no symmetries) and K{t) — 2r 
for all chaotic (n — 0) cavities @. Finally for t > rw 
{th — Pv{v) is the Heisenberg time), K{t) — th |6|. 
Since our purpose is only to estimate the magnitude of 
the fluctuations we follow 9] and assume that the non 
universal part of K{t) is fully described by the contribu- 
tion of the shortest periodic orbit. Within this approxi- 
mation it is straightforward to show that, 



(i?2 



''nun 



L6 



■log(T^f/r„ 



(7) 



We note this result is valid provided that Zgsc 3> Imin- 
Otherwise one would have to multiply the contribution 
of each periodic orbit by its probability to stay inside the 
cavity. 

A similar calculation of the fluctuations 5u{v) of the 
density of energy u{i>) with respect to u{v{v) in a cavity 
with no symmetries shows that, 

(5u(j/))2 _ c2 



""esc 



The dependence of the fluctuations with the cutoff dis- 
tance lesc is expected since the sum over orbits in m(i^) 
Eq. ([31), unlike R{T) Eq.®, does not really converge. 
Without a cutoff u{v) is a series of isolated peaks at the 
natural frequencies of the cavity. 



AN EXAMPLE: A SPHERICAL BLACKBODY. 

The case of a spherical blackbody {n = 3) is of spe- 
cial interest since it can be solved exactly. In addition 
finite size corrections are stronger due to the high degree 
of symmetry of the sphere (the study of a cubic black- 
body was carried out in Ref . [10|] , the non oscillatory cor- 
rections in other simple geometries were investigated in 
ll|). Small multipolar corrections to the spherical shape 
12( 1 are described with the same periodic orbits that in 
the spherical case but adding an additional cutoff in term 
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FIG. 1: Density of energy u(u) in J jrrC' of a blackbody versus 
frequency. The solid line is the Planck's law. The dashed line 
is the analytical prediction for a spherical blackbody Eq.((3]) 
and Eq. (|9]). Finite size effects are clearly visible. The pa- 
rameters chosen (ro = 2cm, T = 5if) fall within the region 
accessible to experiments (see text). In order to reproduce FI- 
RAS resolution only rays shorter than Icut — 2 x W^c/u are 
included in Eq. Q. In the inset we plot Rsca{x) = R{T)/aT'^ 
with X = TrokB /he. 

For this geometry the density p(i^) [4] is known ex- 
plicitly. The closed periodic orbits are given by planar 
regular polygons with a even number of vertexes along 
a plane containing the diameter. No polarization angle 
appears in our expressions due to the planarity of the 
orbits. The length Lp t of the trajectories is given by 
Lp^t = 2pro sin(0) where vq is the radius of the sphere, p 
is the number of vertexes of the polygon and 4> = T^t/p 
with t being the number of turns around the origin of a 
specific periodic orbit. With these definitions, 



p{iy) = pviiy) + PC + K^') 



with 



V 



1^ sm{8nvtro / c) 



c ^ crotTT 



(8) 



(9) 



67r 



t,p>2t 



3/2 



/sin(0) _^ 3/2)7r/2 -I- 2TTvpro sin(0)/c) 

y roTT'^p 

and PC = The sum above will be effectively cutoff 
due to the finite aperture that make the cavity open. 
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We postpone its study to the next section in which we 
address the experimental detection of these effects. 

Corrections to the Stefan-Boltzmann's law can be eval- 
uated explicitly by integration over v. In the limit 

7 flC ^ 7 



R(T) = aT^ 



where &i ~ (due to the fact that only orbits with 
an even number of vertexes contribute to the sum), 
6o ~ sil^T St ~ 2 X 10^** . The second term cor- 
responds to the curvature contribution pc. The last two 
terms come from the oscillating part of the spectral den- 
sity p{v). The sum is dominated by the first terms cor- 
responding to the contribution of the shortest periodic 
orbits and consequently it is quite insensitive to the cut- 
off le.sc ^tnin- 

In order to examine the importance of these finite size 
effects we rescale R{T), 



Rsca{x) 



RjT) 



= 1 



606n 



(11) 



with X = r^TkB / ch. This is an expansion in \/x valid 
for a; ^ 1. The most interesting region (see inset Fig 1.) 
to detect finite size corrections is thus x ^ 1. We note 
that since we go up to fourth order the expansion is still 
accurate for x ^ 2. Values of rg, T compatible with this 
condition are within but close to the limit accessible to 
experiments (ro > 0.5cm and T > 2K). We note that 
a) The dependence of R{T) with T and L for a spher- 
ical cavity Eq. (llOp agrees with the general results ob- 
tained earlier Eq. (O . Only the dimensionless coefficients 
in front of Eq. ([5]) are cavity dependent, b) a similar ex- 
pansion of Rsca valid for cc < 1 can also be carried out. 
However we could not find a region of parameters (T, tq 
with ro > c/v) accessible to experimental verification. 



EXPERIMENTAL DETECTION. 

We determine the range of parameters in which is feasi- 
ble to build a spherical blackbody and examine whether 
finite size corrections in this region could be detected 
with the current technical capabilities. 

We first ask about the optimal radius for experimen- 
tal verification. In general, the larger the radius the more 
difficult it is to keep constant and stable the temperature. 
This constraint excludes cavities with radius tq ^ 0.2m. 
On the other hand in cavities with 1cm it is hard 

to control the geometrical shape. Furthermore, the flux 
of energy through the aperture is too small to be mea- 
sured with the accuracy needed. The optimal radius can 
be estimated to be ro — 0.5 — 7cm. We look at wave- 
lengths X = c/v {ro ^ 8 — 30A) such that the semiclas- 
sical formalism is accurate, namely, — <C 1. Finally 



the optimal temperature is such that for the range of 
wavelengths of interest, A ^ 0.5 — 2mm, the density of 
energy m(^) is a maximum. These conditions are satis- 
fied for T ~ 2.5 — 6K which is accessible by cooling the 
cavity with Helium liquid. Thus we propose that the op- 
timal experimental setting is a spherical blackbody with 
ro - 0.5-7cm,T< 5K, and = c/A ~ 1.5-6x 10"s-^ 
In addition in order to observe the effect of the periodic 
orbits {lesc ^ Imin) the area of the aperture must verify 

In order to proceed we must choose an apparatus of 
measurement. For our purpose it is of paramount im- 
portance that the radiation of the blackbody at T ~ 5K 
can be measured with the highest precision. We also re- 
quest that the resulting flux of energy Eq. ^ through the 
aperture can be measured with a precision much larger 
than the strength of the predicted finite size corrections. 
Finally we require the apparatus to be precise enough 
to discern comparatively close wavelength. In order to 
observe fluctuations in u^v) it is required that the appa- 
ratus can discern frequencies A (A) comparable with the 
contribution to u(iy) of the shortest periodic orbit. In our 
case this corresponds to A(A)/A < 0.05. 

After a careful research, we conclude that instruments 
designed to measure the cosmic microwave background 
are the best suited for our purpose. We focus our analy- 
sis on FIRAS apparatus on board of COBE satellite 
though other differential devices to measure the cosmic 
microwave background could be used as well. By dif- 
ferential it is meant that an internal calibrator nulls the 
external signal coming either from the sky or from the ex- 
ternal calibrator. This mechanism is responsible for the 
FIRAS ability to detect small deviations from a black- 
body source (ICAL in FIRAS). The FIRAS spectrome- 
ter is designed to measure deviations from a blackbody 
spectrum in the region A = 0.1 — 10mm. The intrinsic 
frequency resolution is A(i/)/i/ ^ 5 x lO^'^. This im- 
plies that the contribution of periodic orbits (light rays) 
Ip ^ Icut = 2x lO^c/i^ will not be detected by the appara- 
tus. In the range of frequencies of interest the equivalent 
noise power 4 x 10~^^W/^/ Hz is less than 1% of the 
measurement. The blackbody temperature can be set in 
the range T — 2 ~ lOK with a precision of a few rriK. 
In Fig.l we plot u{i') for a spherical cavity for a set of 
parameters (ro = 2cm, T = 5K and Icut = 2 x lO^c/iy) 
accessible to experimental verification. From the figure 
it is clear that an apparatus with FIRAS specifications is 
capable to detect finite size corrections in spherical black- 
body cavities. Similar results will be obtained for other 
symmetric cavities such as cylinders or cubes. 

We note that a condition for the experimental obser- 
vation of these effects is that the radiation is coupled out 
of the cavity in such a way the original modes of the cav- 
ity are not seriously affected. Previously we have shown 
that the aperture of the blackbody acts just as a natural 
cutoff for long periodic orbits. Therefore this is not a 
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problem for the experimental verification of our results. 
Indeed in recent experiments 0| it has been possible to 
measure the natural modes of a cavity with a microwave 
source with a precision enough to even test quantitatively 
semiclassical estimation for the number of modes. 



APPLICATIONS: SONOLUMINESCENCE AND 
THE COSMIC MICROWAVE BACKGROUND. 

Differential apparatus as FIRAS compare the sky sig- 
nal with some calibrator on board which is considered a 
'perfect' blackbody. These calibrators are effective black- 
body cavities which may be affected by the finite volume 
corrections reported in this letter. These corrections may 
very well be confused with those coming from the sky. 
Thus on board calibrators (or cavities used as reference 
sources in the measurement, for instance, of the spec- 
tral emissivity of solid samples should be carefully 
designed to minimize these effects. 

We speculate that finite size corrections may also play 
a role in cosmological problems. The CMB observed to- 
day had its origin in photons from the last scattering 
surface. For angular separations larger than 1.5 degree, 
temperature fiuctuations in the CMB comes from regions 
that have never been in thermal contact. The origin of 
these fluctuations is supposed to date back to the infla- 
tion time. We speculate that flnite size effects may play 
some role in the generation of these fluctuations. In order 
to determine the true relevance of this effect it would be 
necessary to: a)understand the physical mechanisms such 
as decoherence that may lead to a suppression of these ef- 
fects, b) study these deviations in a non Euclidean cavity 
resembling the global geometry of the universe. It is en- 
couraging that in Euclidean cavities it can be shown that 
the shape of the finite size corrections to the blackbody's 
laws are not modified as universe expands. 

Sonolumincsccnce, the transduction of sound into 
light, [isl. [l6|. occurs when the pulsations of an almost 
spherical bubble of gas in water (or other substances [14]) 
produced by a standing sound wave attains sufficient am- 
plitude so as to emit periodic picosecond flashes of light. 
In a certain range of frequencies the spectrum is consis- 
tent with a blackbody 14, of typical size tq ~ 1/im 
and temperatures T ~ 10*K (for bubbles of Xenon in 
sulphuric acid rp ~ 3.9/im and T ^ TOOOi^T). In this 
range of parameters finite size effects should be observ- 
able. However we note in this case the blackbody is not 
the perfecting conducting cavity with an aperture stud- 
ied previously but rather a material of refraction index ni 
immersed in a medium with refraction index n2. A nat- 
ural question to ask is whether the analytical results of 



the previous section are also applicable in this situation. 
The answer is affirmative. The only difference is the way 
in which the cutoff of long periodic orbits is defined. For 
TT-i > n2 the blackbody is effectively an open cavity. A 
given periodic orbit (light rays) inside the material of re- 
fraction index ni will contribute to the spectral density 
provided that the angle of incidence 6 > dc (measured 
from the normal) with 0c = arcsin(n2/ni). Periodic or- 
bits such that 6 < dc do not contribute since the light 
ray escapes to the medium with n2 (l7l |. 

In conclusion, by using semiclassical techniques, we 
have derived explicit expressions for the flnite size cor- 
rections of the blackbody radiation's laws as a function 
of the temperature, size and shape of the cavity and area 
of the aperture. We have also shown that the experimen- 
tal detection of these finite size corrections is within the 
reach of current experimental capabilities and may be of 
relevance in sonoluminescence. 
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acknowledges financial support from a Marie Curie Out- 
going Action, contract MOIF-CT-2005-007300. 



[1] 
[2] 
[3] 
[4] 
[5] 

[6] 

[7] 
[8] 
[9] 

[10] 
[11] 



[12] 

[13] 
[14] 
[15] 



[16] 
[17] 
[18] 



R.Balian and B. Duplantier, Ann. Phys. 104, 300 (1977). 
R. Balian and C. Bloch, Ann. Phys. 64, 271 (1971). 
R. Balian and C. Bloch, Ann. Phys. 60, 401 (1970). 
R. Balian and C. Bloch, Ann. Phys. 69, 76 (1971). 
M. Brack, R.K. Bhaduri, Semiclassical Physics, Addison- 
Wesley, New York, (1997). 

P. Leboeuf, Lect. Notes Phys. 652, Springer, Berlin Hei- 
delberg 2005, p. 245, J. M. Arias and M. Lozano (Eds.). 
C. Dembowski, et al., Phys. Rev. Lett. 89, 064101 (2002). 
R.A. Jalabert, cond-mat/99 12038 

H. 01ofsson,et al., arXiv:0704.2310 ,0. Bohigas and P. 
Leboeuf, Phys. Rev. Lett. 88 (2002) 092502. 
K.M. Case and S.C. Chin, Phys. Rev. A 1 1170 (1970);H.- 
T. Elze and W. Greiner, Phys. Lett. B179 (1986) 385. 
H.P. Baltes and F.K. Kneubuhl, Phys. LettSOA, 360 
(1969);H.P. Baltes and F. K. Kneubuhl, Hel. Phys. Acta 
45, 481 1972. 

P. Meier,et al., Z. Phys. D 41 281 (1997). 

J.C. Mather,et al., Ap.J 512 511 (1999). 

S. D. Hopkins et al., Phys. Rev. Lett. 95, 254301 (2005). 

G. Vazquez,et al., Phys. Rev. Lett. 88 197402 (2002);G. 

Vazquez, C. Camara, S. Putterman and K. Weninger, 

Opt. Lett. 26 575 (2001). 

B.P. Barber, et al., Phys. Rep. 281 66 (1997). 

M. Lebental, et al., Phys. Rev. A 76, 023830 (2007). 

N. Pavloffand C. Schmit, Phys. Rev. B 58, 4942 (1998). 



